Theory of RF-Spectroscopy of strongly interacting Fermions 
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We show that strong pairing correlations in Fermi gases lead to the appearance of a gap-like 
structure in the RF-spectrum, both in the balanced superfluid and in the normal phase above the 
Clogston-Chandrasekhar limit. The average RF-shift of a unitary gas is proportional to the ratio 
of the Fermi velocity and the scattering length with the final state. In the strongly imbalanced 
case, the RF-spectrum measures the binding energy of a minority atom to the Fermi sea of majority 
atoms. Our results provide a qualitative understanding of recent experiments by Schunck et.al. 

PACS numbers: 03.75.Ss, 05.30.Fk, 32.30.Bv 



According to the fifty year old microscopic theory of 
Bardeen, Cooper and Schrieffer, the phenomenon of su- 
perfluidity in a system of fermions is connected with 
the formation of bound pairs. In the weak coupling 
limit, where the formation of pairs and their condensa- 
tion appears simultaneously, the transition to the super- 
fluid state is associated with the appearance of a gap in 
the fermionic excitation spectrum. For strong coupling, 
however, this simple connection is no longer valid and 
bound pairs of fermions may exist even in the normal 
state. This phenomenon is well known from the pseudo- 
gap phase in high temperature superconductors, where 
a d-wave pairing gap appears on the Fermi surface at 
temperatures far above the superconducting transition 
temperature A much simpler example is realized by 
ultracold fermions near a Feshbach resonance, which pro- 
vide a perfectly controllable model system to study the 
effects of strong pairing interactions ^] . In the case of an 
equal population of the two hyperfine states undergoing 
pairing, the ground state is superfluid at arbitrary val- 
ues of the scattering length. A microscopic signature of 
pairing in ultracold Fermi gases has first been obtained 
by Chin et.al. -S] through RF-spectroscopy. The RF- 
field drives transitions between one of the hyperfine states 
1 2) = U) which is involved in the pairing and an empty 
hyperfine state |3) which lies above it by an energy hLU23 
due to the magnetic field splitting of the bare atom hy- 
perfine levels. In the absence of any interactions, the 
spectrum exhibits a sharp peak at w = W23. Pairing be- 
tween the two lowest hyperfine states |1) and |2) leads to 
an upward shift of this resonance. The shift essentially 
follows the two-particle binding energy on the BEC-side 
of the crossover but stays finite on the BCS-side, where 
the appearance of a bound Cooper pairs is a many-body 
effect [3| . A theoretical explanation of these observations 
can be given by extending the BCS description of pairing 
to the strong coupling regime and neglecting interactions 
involving state |3) [J HI. In a homogeneous system, the 
resulting RF-spectrum exhibits a peak at energies around 
A^/ /i, which is of the order of the energy gap A w 0.5 ep 
at the unitarity point. Since pairing appears already in 
the normal state above Tc, the RF-shift does not di- 



rectly measure the superfluid order, however The 
importance of understanding the relation between RF- 
spectra and the nature of the many-body states involved, 
is underlined by recent experiments in imbalanced gases 
Q. There, a shift in the RF-spectrum is observed which 
hardly changes between the balanced superfluid and a 
normal ground state beyond a critical population imbal- 
ance, where superfluidity is destroyed by a sufficiently 
large mismatch of the Fermi energies even at T = (this 
is the analog of the Clogston-Chandrasekhar limit in su- 
perconductors). In this work, we present a theory of 
RF-shifts in both balanced and imbalanced Fermi gases, 
which provides a qualitative understanding of these ob- 
servations. In particular, we show that the average fre- 
quency shift in the balanced superfluid at unitarity (i.e. 
at inflnite scattering length) is linear in the Fermi ve- 
locity and inversely proportional to the scattering length 
ai3. In the non-superfluid state beyond the Clogston- 
Chandrasekhar limit, pair fluctuations give rise to sharp 
peaks in the RF-spectrum which are associated with the 
binding of t |-pairs even in the absence of long range 
phase coherence. 

Within linear response theory, which is adequate for 
RF pulses short compared to the Rabi oscillation period 
of the bare 2-3 transition, the number of particles trans- 
ferred from state |2) to state |3) per unit time is given 
by 
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where lo = ujl — W23 denotes the detuning of the RF field 
from the bare 2-3 transition. Since particles in state |3) 
have a nonvanishing interaction with those in states |1) 
and |2) Q, the response function in equation (1) does 
not factorize into one particle functions, making a full 
calculation of the spectrum very difficult. Nevertheless, 
near T — where only a single peak is observed in the 
RF-spectrum, its position can be determined from a sum 
rule approach Q. In particular, the first moment uj = 



2 



J dLuujI{uj)/ J dujl(uj) is given by 
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Here i/^s ^-nd denote the interaction Hamihonians 
between the respective states, while N2 and N3 denote 
the total number of particles in states |2) and |3). The 
gtj are the bare interaction constants arising in the pseu- 
dopotential interaction Hamiltonian 



g,, / d^x4{x)i;]{x)iPj{x)A{x). (3) 



They are related to their renormalized values gij 
ATrfi^aij /m by 
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where the s-wave scattering lengths between states 

i and j, m is the mass of the particles and Ek = fi^k^ /2'm 
the single particle energy. Note that the interaction 523 
between states 2 and 3 drops out quite generally, because 
and Hyif commute. Moreover, there is no shift of the 
RF peak if the interaction strengths 912 and (713 are equal, 
a case, where all interaction effects are cancelled exactly 
^,9]. Since (H-^^^) is of order A^3, the first term in ^ is 
negligible compared to the second term if N2 ^ N3. The 
average shift of the RF-spectrum then simplifies to 
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Here, the second form is obtained by expanding 1 — 
913/912 to leading order in the upper cutoff A of the 
momentum integral in (4). Evidently, for vanishing in- 
teractions gi3 = gi3 = with state 3, the RF-shift just 
measures the (negative) interaction energy per particle 
in the state 2. Within a pseudopotential description, 
however, the interaction energy {H12) ~ A diverges lin- 
early with the cutoff. It is thus sensitive to the range 
of the interactions, which is set equal to zero in the 
pseudopotential. In terms of the spectrum I{lo), this 
divergence shows up as a slow decay ^ uj~^/^ at 

large frequencies, leading to a divergent first moment, as 
is easily seen within a BCS-description with a constant 
gap A. Remarkably, for finite interactions gis 7^ 0, the 
second form of (5) gives a result for the frequency shift 
which is well defined and finite in the limit A ^ 00. As 
shown by Tan the total energy of the balanced gas 
can be obtained from the momentum distribution via 
E = 2^j,ek("'k — C/k'^) up to a constant, which is ir- 
relevant for the calculation of the limit Here 
C is the constant arising in the asymptotic behavior 
limrik — C/k'^ of the momentum distribution at large 
momenta. Evidently, the interaction contribution to the 



total energy is just (^12) = ~2C^^. ek/fc"* — CA. In 
troducing a dimensionless constant s via C 
shift of the RF-spectrum 



sk'^, the 
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of the balanced gas is completely determined by the uni- 
versal constant s, the Fermi energy ep = h'^kp/{2m) of 
the balanced, non-interacting gas and the renormalized 
interaction constants 312 and gi3. The expression is fi- 
nite for all coupling strengths gi2 and evolves smoothly 
from the BCS- to the BEC-limit. Within an extended 
BCS-description of the ground state wavefunction, the 
product s'"' • 4e|, = is precisely the square of the 
gap parameter. In weak coupling, our result then coin- 
cides with that obtained by Yu and Baym [§| , except for 
the mean field shift, which is not contained in the re- 
duced BCS Hamiltonian. In the BEC-limit, where the 
BCS-groundstate becomes exact, the asymptotic behav- 
ior Abec = 4£f/ VSTrftpo^ gives fuv = 2£b(l - 012/013), 



where et — h^/r 



^12 



is the two-particle binding en- 



ergy. It is straightforward to show, that this is pre- 
cisely the average shift for bound-free transitions fol- 
lowing from a detailed calculation of the RF-spectrum 



in the molecular limit by Chin and Julienne The 
most interesting regime is that around the unitarity limit 
1/512 = 0. At this point, the average RF-shift is given 
by Q = —OAQvp/ais, which varies like the square root of 
the Fermi energy ep = mvp/2. The constant s = 0.098 
is obtained from the recent calculations of the crossover 
thermodynamics by Haussmann et al. [13] ■ Our result 
for the homogeneous gas can be compared directly with 
locally resolved RF-spectra by Shin et al. [3] ■ Account- 
ing for the enhancement of the local Fermi velocity at 
the trap center by a factor « 1.25 due to the attractive 
interactions, the predicted average shift a; = 27r- 28.9kHz 
[3| is considerably larger than the measured position of 
the peak near 15 kHz. This is probably due to the fact, 
that Lu has a considerable contribution from the higher 
frequency part of the spectrum. A crucial prediction of 
our theory is the linear behaviour of the average RF-shift 
with the Fermi momentum. Experimentally, the spatial 
resolution necessary to distinguish this from the naive 
e_F-scaling has not yet been achieved [l3|- 

To discuss the situation with a finite imbalance, it is 
convenient to introduce two distinct chemical potentials 
for the states undergoing pairing, defined hy = jj, + h 
and /i| ~ 11 — h. Since the ground state of the spin 
balanced gas is a superfluid with a gap for fermionic ex- 
citations, it will be stable over a finite range h < he 
of the chemical potential difference. In the BCS limit, 
the associated Clogston-Chandrasekhar critical field he = 
Abcs/v^ is exponentially small. Near the unitarity 
point, the absence of a second energy scale implies that 
the critical field he beyond which a non-zero polariza- 
tion appears, is on the order of the bare Fermi energy 





FIG. 1: Lowest order self energy diagram for the minority 
component (|) Green's function and Bethe-Salpeter equation 
for the vertex part in ladder approximation. 



Sf of the balanced two-component Fermi gas. From 
fixed node diffusion Monte Carlo calculations the result- 
ing numerical value in the continuum case at unitarity 
he — 0.96 /i « OAep 15|. The phase for h > he is 



IS 

a non-superfluid, polarized mixture of the different spin 
states. For large enough fields, the system will eventually 
be completely spin polarized. At unitarity, the associated 
saturation field hg was determined by Chevy [16] using 
a variational calculation of the energy change asso- 
ciated with adding a single |-particle to a Fermi sea of 
"f-particles. This leads to an upper bound fi^ < —0.60 /i| 
at the unitarity point, where = 2^^^eF is the Fermi en- 
ergy of the completely spin polarized gas. The saturation 
field thus obeys the inequality hg > 0.8 /i| — 1.27 ep- At 
unitarity, therefore, there is a wide regime he < h < hs 
of an intermediate phase between the balanced super- 
fluid and a fully polarized gas. While superfluidity is 
quenched in this phase, the strong interactions between 
particles in states |1) and |2) still give rise to large fre- 
quency shifts in the RF-spectrum, as will be shown be- 
low. To study the effect of pairing fluctuations on the 
imbalanced Fermi gas above the Clogston-Chandrasekhar 
limit, we calculate the pair-fluctuation spectrum from the 
two-fermion Green function, using a non-selfconsistent T- 
matrix app roach, similar to the approach by Combescot 



et al. 17|. Such a perturbative analysis is reasonable. 



since the states which are coupled through the interac- 
tion Hamiltonian are separated by an energy gap of width 
2h. A usual ladder approximation is used to incorporate 
the effects of the attractive t i-interaction on the vertex 
part, whereas the self energy is calculated at the one- loop 
level, including vertex corrections. The basic equations 
for the polarization loop L, vertex part F and self-energy 
part of the minority species are given by (we take 
units in which h = 1) 



iiln - iUJ,n) 



1 



l/.g-i(q, i^r^ 



(7) 

(8) 




FIG. 2: Pair binding energy in units of /if at unitarity 
as a function of h at T=0. For h > ji the binding energy is 
constant and given by 0+ = 0.61/it. 
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where Q^^^ and Cjj"'' are the bare Matsubara-Green's 
functions of the majority- and minority component and 
(3 = l/kgT is the inverse temperature. aj„ = {2n + l)T:/ (3 
and ri„ — 2Tin/j3 with n & denote fermionic and 
bosonic Matsubara frequencies respectively. After eval- 
uating the Matsubara summation and analytic continu- 
ation, the vertex part can be calculated analytically at 
T = 0. In the regime h > fi (i.e. essentially beyond 
the Clogston-Chandrasekhar field h > he ^ 0.96/x), one 
obtains for q = 0, a; > — 2/i 
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where 0{x) is the unit step function and kp^- is defined 
via kp]' — ^2m[i\ jh. For h > fi the retarded ver- 
tex r^(q = Ojw) has a single pole on the real axis at 
uJq — 2h — r2+ with $7+ > (note that for h < ^ the 
vertex has two real poles) . Physically, this pole describes 
an excitation in which two fermions with opposite spin 
and vanishing total momentum form a pair at the Fermi 
energy of the majority component with binding energy 
57+ . A similar structure was first discussed for weak cou- 
pling by Aleiner and Altshuler [13] in the context of small 
superconducting grains. Remarkably, as shown in Fig. 
[21 the pair binding energy in units of /i| is constant for 
h > fj. and agrees well with the value 0.6/i| for the binding 
energy of a single down spin in the presence of a Fermi 
sea of majority atoms as calculated by Chevy [16| . The 
retarded self energy for the minority component in the 
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FIG. 3: RF spectra at unitarity for different imbalances at 
T = (Intensity in arbitrary units; he = 0.96/i). 



normal state is given by 

d^q dz 



X Imr^(q, z) - npiz) ImG2_^^(q - k 



X r^(q,z + w)| 



(11) 



with ns and np denoting the Bose- and Fermi- 
distributions. This result enables us to calculate RF spec- 
tra explicitly in the limit of vanishing (723 and (713, where 
the expactation value in equ. Ilj can be factorized. In 
this case, one obtains 0, 0| 



d^k 

^ImGf(k,ek 



if state 1 3) is initially empty. In figure [3] we have numer- 
ically evaluated the resulting RF-spectra at unitarity for 
different fields above the Clogston-Chandrasekhar limit. 
The calculation explains two features which are seen in 
the experimental data namely the shift of the RF 
peak due to pairing fluctuations in the normal state and 
the decreasing linewidth with increasing population im- 
balance. The onset of the RF-spectrum coincides with 
the pair binding energy ri_|_ « 0.6 /Li| for h > he, which is 
independent of the imbalance. In the presence of a finite 
|1) — 1 3) interaction, the detailed spectrum I{uj) can not 
be calculated analytically. Its first moment, however, is 
again determined by the sumrule equ. ([5]). Evaluating 
the interaction energy (iJj^j) using the variational wave- 
function of Chevy [l6|, it turns out that the resulting 
average RF-shift for an almost completely polarized gas 
is equal to w = —O.iAhkp^ /mais. Due to the sharp- 
ness of the peak in this limit, the average shift in the 
strongly imbalanced gas coincides with the experimen- 
tally observed peak position. For the parameters in [l3| , 
we obtain an average RF-shift (D = 27r • 17 kHz at the 



trap center for strong imbalance, close to the observed 
value in the balanced case. Our theory thus accounts for 
the observation by Schunck et al. Q, where an average 
over the trap is involved, that there is hardly any differ- 
ence in the RF-shift between the balanced and strongly 
imbalanced gas. 

In conclusion, we have given a theory of RF-spectra 
in ultracold Fermi gases which includes interactions be- 
tween all three states involved. In the balanced unitary 
gas, the average RF-shift is proportional to —svp/ais, 
where s is a universal constant characterizing the fermion 
momentum distribution at large wave vectors. In the im- 
balanced case, the RF-spectrum exhibits a sharp peak 
arising from the binding energy of a t i-pair which is fi- 
nite even in the non-superfluid state. Including a finite 
value of ai3, the resulting average shift is close to the 
peak shift in the balanced case. 
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with W. Rantner, Yong-il Shin and M. Zwierlein. This 
work was supported by the DFG Forschergruppe " Strong 
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Note added in proof: Equivalent results for the RF- 
shift of balanced gases have been obtained independently 
by Baym et al. [l9|. In fact our value for the prefactor 
'm Cu — —OAGvp/ais agrees well with the value obtained 
in this reference, using a different method. 
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